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Linking two consecutive non-merging magnetic clouds with their 
solar sources 

S. Dasso^'^, C.H. Mandrini^, B. Schmieder^, H. Cremades^, C. Cid^, Y. Cerrato^, 
E. Saiz^, P. Demoulin^, A.N. Zhukov^'^, L. Rodriguez^, A. Aran^, M. Menvielle^, 
and S. Poedts^° 

Abstract. On May 15, 2005, a huge interplanetary coronal mass ejection (ICME) was 
observed near Earth. It triggered one of the most intense geomagnetic storms of solar 
cycle 23 {Dstpeak = -263 nT). This structure has been associated with the two-ribbon 
flare, filament eruption and CME originating in active region (AR) 10759 (NOAA num- 
ber). We analyze here the sequence of events, from solar wind measurements (at 1 AU) 
and back to the Sun, to understand the origin and evolution of this geoeffective ICME. 
Prom a detailed observational study of in situ magnetic field observations and plasma 
parameters in the interplanetary (IP) medium, and the use of appropriate models, we 
propose an alternative interpretation of the IP observations, different to those discussed 
in previous studies. In our view, the IP striicture is formed by two extremely close con- 
secutive magnetic clouds (MCs) that preserve their identity during their propagation through 
the interplanetary medium. Consequently, we identify two solar events in Ha and EUV 
which occurred in the source region of the MCs. The timing between solar and IP events, 
as well as the orientation of the MC axes and their associated solar arcades are in good 
agreement. Additionally, interplanetary radio type II observations allow the tracking of 
the multiple structure through inner heliosphere and to pin down the interaction region 
to be located midway between the Sun and the Earth. The chain of observations from 
the photosphere to interplanetary space is in agreement with this scenario. Our anal- 
ysis allows the detection of the solar sources of the transients and explains the extremely 
fast changes of the solar wind due to the transport of two attached — though non-merging — 
MCs which affect the magnetosphere. 



1. Introduction 

Coronal mass ejections (CMEs) remove plasma and mag- 
netic field from the Sun expelling them into interplanetary 
space. When they are detected in the interplanetary (IP) 
medium, they are called interplanetary coronal mass ejec- 
tions (ICMEs). The transit time of an ICME from the Sun 
to 1 AU is typically in the interval 1 to 5 days [Gopalswamy 
et al, 2000,' 2001a; Rust et al, 2005]. A subset of ICMEs, 
called magnetic clouds (MCs), arc characterized by an en- 
hanced magnetic field strength, a smooth and large rotation 
of the magnetic field vector, and low proton temperature 
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[Burlaga et al, 1981]. Fast and large magnetic clouds are 
mostly observed at 1 AU in the declining phase of a solar 
cycle, as it happened in the last one during 2003-2005 [e.g., 
Culhane and Siscoe, 2007]. 

Interplanetary type II bursts permit the tracking of fast 
ICMEs through the analysis of radio emission frequency drift 
[Reiner et al., 1998, 2007; Hoang et al., 2007]. Assuming the 
type II emission to be produced at the fundamental or sec- 
ond harmonic of the local plasma frequency, and with help 
of a heliospheric density model, the radio frequency can be 
converted to radial distance. This method has still limita- 
tions and allows for possible different interpretations due to 
the patchiness and the frequency range of the observed radio 
emissions. The evolution of an ICME can also be followed 
along its journey using a combination of solar and IP obser- 
vations together with the results of numerical simulations 
[e.g., Wu et al, 1999]. 

Coronal mass ejections are frequently associated with fil- 
ament eruptions. The directions of the MC axes are found 
to be roughly aligned with the disappearing filaments [Both- 
mer and Schwenn, 1994, 1998], preserving their hehcity sign. 
This result has been found for some individual cases by 
Marubashi [1997]; Yurchyshyn et al. [2001]; Ruzmaikin et al. 
[2003]; Yurchyshyn et al. [2005]; Rodriguez L. et al. [2008]. 
Some quantitative (quantifying magnetic fluxes and holici- 
ties) studies of MCs and their solar sources have been also 
done [Mandnm et al., 2005; Luoni et al., 2005; Longcope 
et al., 2007]. But again, it is not so easy to quantify this 
association and further developments are needed to really 
understand the fundamental transport mechanisms and in- 
teraction with the ambient solar wind (see, e.g., the review 
by DemouUn [2008]). 

Furthermore, when multiple CMEs are expelled from the 
Sun, they can be merged leading to the so-called "CME 
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cannibalism" [Gopalswarny et al., 2001b], mainly as a conse- 
quence of magnetic reconnection. The interaction of two 
CMEs in favorable conditions for reconnection has been 
studied by Wang et al. [2005]. However, from a theoretical 
point of view, two original structures can be preserved with 
an appropriate orientation of the ejected flux rope (yielding 
almost parallel interacting magnetic fields) . For example, in 
the MHD simulations of Xiong et al. [2007] the two inter- 
acting flux ropes preserved their identity while evolving in 
the IP medium until they reach the Earth environment or 
even beyond. 

In the magnctohydrodynamic (MHD) framework, a mag- 
netic cloud configuration in equilibrium can be obtained 
from the balance between the magnetic Lorentz force and 
the plasma pressure gradient. Several magnetostatic models 
have been used to describe the configuration of MCs. Fre- 
quently, the magnetic field of MCs has been modeled by the 
so-called Lundquist's model [Lundquist, 1950], which consid- 
ers a static and axially symmetric (cylindrical) linear force- 
free magnetic configuration neglecting the plasma pressure 
[e.g., Goldstein, 1983; Burlaga, 1988; Lepping et al, 1990; 
Lynch et al, 2003; Dasso et al., 2005; Leitner et al., 2007; 
Xiong et al., 2007]. The azimuthal and axial field compo- 
nents of the flux rope in this classical configuration are de- 
fined by — B[)Ji{cvr) and Bz = BoJo{cir), where r is 
the distance to the MC axis, J„ is the Bessel function of 
the first kind of order n, Bq is the strength of the field 
at the cloud axis, and a is a constant associated with the 
twist of the magnetic field lines. Some other refined models 
have also been proposed to describe the magnetic structure 
of clouds [e.g., Hu and Sonnerup, 2001; Vandas and Ro- 
mashets, 2002; Cid et al., 2002]; in particular, some of them 
consider an elliptical shape [Hidalgo et al., 2002a; Hidalgo, 
2003] that allows the description of possible distortions of 
the structures. 

On May 15-17, 2005, the strongly southward interplan- 
etary field (above 40 nT) and the high solar wind velocity 
(close to 1000 km s^^), observed by the Advanced Com- 
position Explorer (ACE), arc the cause of a super geomag- 
netic storm with a depression of the Dst index reaching 
-263 nT. The IP structure has the characteristics of a MC. 
Yurchyshyn et al. [2006] associated this structure with the 
two-ribbon M8.0 X-ray class flare on May 13 at 16:32 UT in 
AR 10759, accompanied by a filament eruption and CME. 
These authors identified a single MC, starting at ~ 00:00 UT 
on May 15 and ending at ~ 09:00 UT on May 16. An 
event with different boundaries, much smaller than the one 
described by Yurchyshyn et al., was identified by R. Lop- 
ping (start on May 15 at 05:42 UT and end on May 15 at 
22:18 UT). The event was qualified with intermediate fitting 
quality (quality 2) and had a large impact parameter (cata- 
log at http://lepmfi.gsfc.nasa.gov/mfi/mag_cloud_Sl.html). 

The solar activity evolution along May 13 has been well 
described by Yurchyshyn et al. [2006] and, more recently, 
by Liu et al. [2007]. Both papers refer mainly to the M8.0 
event and conclude that the eruption of a large sigmoidal 
structure launches the CME observed by the Large Angle 
and Spectroscopic Coronagraph [LASCO, Brueckner G. E. 
et al, 1995] at 17:22 UT. 

In this paper, we present a detailed analysis of the solar 
wind conditions at the Lagrangian LI point using in situ 
observations (magnetic field and bulk plasma properties) 
and propose an alternative interpretation of the IP obser- 
vations. In our view, the observed ICME is in fact formed 
by two MCs. Since a single solar event cannot explain the 
arrival of two MCs at 1 AU, we revisit the evolution and 
activity of AR 10759 and other regions present in the solar 
disk. We start searching in the Sun for two possible sources 
of the clouds, even earlier than May 13. We identify a pre- 
vious candidate event in Ha and EUV data same day at 
12:54 UT in AR 10759, which was classified as a CI. 5 flare 



in GOES. These two events, the one at 12:54 UT and the 
one at 16:32 UT, give rise to two two-ribbon flares along 
different portions of the AR magnetic inversion line. The 
orientations of the magnetic fields associated to both solar 
events are in good agreement with the fields observed in 
their associated MCs. 

In Section 2 we analyze the solax wind conditions near 
Earth; in particular, we derive the orientation of the axis 
for the two MCs by fitting a model to the observations. In 
Section 3, we revisit the solar events using ground based 
and satellite data. In Section 4 we present radio type II 
remote observations, which let us track the multiple struc- 
ture through the inner heliosphere. Finally, in Section 5, 
we describe the most plausible physical scenario behind the 
chain of events from May 13 to 17, 2005, and we give our 
conclusions. 

2. A Large Interplanetary Coronal Mass 
Ejection: May 15-17. 

2.1. In situ plasma and magnetic field observations 

The solar wind data were obtained by ACE, located at 
the Lagrangian point LI. The magnetic field observations 
come from the Magnetic Fields Experiment (MAG) [Smith 
et al, 1998] and the plasma data from the Solar Wind Elec- 
tron Proton Alpha Monitor (SWEPAM) [McComas et al, 
1998]. 

Figure 1 shows the magnetic field (in the Geocentric Solar 
Ecliptic, GSE, system) and plasma conditions in the solar 
wind for the long-lasting studied event (00:00 UT on May 
15 to 16:00 UT on May 17, 2005), and the consequent mag- 
netospheric activity {Dst index). 

The data in Figure 1 show several features consistent with 
the presence of ICMEs/MCs: low proton plasma beta, pro- 
ton temperature lower than expected (Tex) for a typical so- 
lar wind at a given observed bulk velocity [Lopez, 1987], an 
almost linear velocity profile, and a smoothly varying mag- 
netic field orientation of high intensity. In particular, from 
vertical dashed line '1' to '10' (for timings of tick numbers 
see next paragraphs and Table 1) the magnetic field inten- 
sity is strongly enhanced, although the scale needed to show 
the whole event in one plot does not allow to observe this 
enhancement at the end of the time interval. The tempo- 
ral length of the event is quite large, so that it corresponds 



Table 1. Timings (dd,hh:mm UT) and substructures inside 
the ICME of May 15-17, 2005,their boundaries are identified 
with numbers in Figure 1. 



Tick number 


Timing 


Substructure 


1 


15,02:11 


sheath 


2 


15,05:42 


MCI 


3 


15,10:20 


backl 


4 


15,14:10 


MC2 


5 


15,21:40 


MC2 


6 


16,00:15 


MC2 


7 

8 


16,04:10 
16,13:00 


back2 
back2 


9 


17,09:37 


back2 


10 


17,10:30 
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Figure 1. In situ (at LI) plasma and magnetic field observations (ACE) of the ICME on May 2005. 
From top to bottom: absolute value of the magnetic field (B = |i?|), magnetic field vector orientation 
(GSE): latitude (Qb) and longitude (0s), bulk velocity (F), expected (continuous line) and observed 
(dots) proton temperature (Tp), proton density (rip), proton plasma beta ((8p), and Dst index. Vertical 
lines mark different interfaces (see Section 2.1 for a description and Table 1 for timings). Horizontal 
dotted lines in Qb, 4>b, Pp, and Dst mark values at 0, 180, 0.1, and -100 as a reference, respectively. 



to a very extended region, ~ 0.93 AU, one of the largest 
ICMEs ever observed (see, e.g.. Table 1 of Liu et al. [2005] 
and Figure 1 of Liu et al. [2006a]). At '10', the magnetic 
field recovers its background value of ~ 5 nT. 

From a comparative study between magnetic clouds and 
complex ejecta, Burlaga et al. [2001] found that while most 
of the clouds were associated with single solar sources, nearly 
all the complex ejecta could have had multiple sources. Evi- 
dence was also found indicating that long duration complex 
merged interaction regions (with radial extents of ~ 0.7 
AUs) can be produced by the interaction of two or more 
CMEs/MCs/shocks [Burlaga et al, 2001, 2003]. 

In front of the ICME, from 00:00 UT to 02:11 UT on 
May 15 (labeled as '1'), the solar wind presents typical con- 
ditions with a value of B 6 nT, an increasing velocity pro- 
file (starting at ~ 400 km s"'^ and reaching ~ 500 km s"'^), 
observed proton temperature similar to that expected one 
{T^x) for a typical solar at same velocity, an increasing pro- 
ton density profile from 3 to 10 protons per cnr\ high values 
of proton plasma /3p (reaching ~ 10-100), and values of the 
Dst index {Dst ~ 0) corresponding to relatively quiet ring 
current conditions. 

Observations at '1' suggest the existence of a strong lead- 
ing edge shock (not fully confirmed from observations of V 



because of a gap in plasma data from 02:11 UT to 03:50 UT 
on May 15). Just behind the shock, from '1' to '2', a typical 
ICME sheath is present (e.g., high level of fluctuations in B, 
enhanced B strength, high mass density, high /3p). 

A structure with a very high B (50-60 nT), one of the 
highest values ever observed in the solar wind at 1 AU, is 
found between '2' and '3'. In this range of time, a large scale 
coherent rotation of the magnetic field vector is present, with 
B going from south to north (see 9b panel). This is a left 
handed SEN flux rope in the classification of Bothmer and 
Schwenn [1998]. Note the presence of a sudden change of 
the sense of rotation together with a magnetic discontinu- 
ity (MD) at '3' (clearly observed in 9b). A current sheet 
(i.e., a discontinuity in the observed time series of the mag- 
netic field vector) is expected to be present at the interface 
that separates two magnetic regions with different magnetic 
connectivity, i.e., with different magnetic stress. Thus, we 
interprete the discontinuity at '3' as the signature of the end 
of a first substructure inside the large ICME. This substruc- 
ture '2-3' has also specific physical properties: there are no 
signatures of expansion (the profile of V does not show a 
significant slope), and the observed temperature decreases 
in time (hotter near the sheath and colder near substructure 
'3-4'). As a consequence of the decreasing Tp profile, values 
of /3p also decrease, being /3p ~ 10~^ near '2' and reaching 
values as low as /3p ~ 10~ near '3'. Thus, the substruc- 
ture '2-3' presents some MC signatures (low Pp, coherent 
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rotation of the magnetic field vector) but not all ICMEs 
features (e.g., Tp is not significantly lower than Te^)- This 
MC presents a short temporal duration with a spatial size 
along the Sun-Earth direction of ~0.09 AU (~4 hs multi- 
plied by ~900 km s~^). We interpret the short duration of 
the first MC as due to the compression made by MC2 at 
its trailing edge; it is expected that the expansion rate (and 
the consequent cooling) of the first cloud diminished during 
part of its journey, where the second MC is pushing it from 
behind [Wang et al., 2003], consistently with observations 

{Tp ~ Tea,). 

After the MD at '3', 6b continues increasing up to '4', 
where there is another (weaker, but significant) MD. The 
region from '3' to '4' shows the characteristics of the back of 
fiux ropes, previously found by Dasso et al. [2006] in a sig- 
nificantly large non-expanding magnetic cloud (Oct. 1995) 
and by Dasso et al. [2007] in another huge magnetic cloud 
(Nov. 2004) in strong expansion. The formation of this back 
feature is a consequence of previous magnetic reconnection 
between the front of the flux rope and its environment; mag- 
netic fiux is removed from the flux rope front (the front is 
peeled) while its counterpart in the rear part still remains 
(Figure 6 of Dasso et al. [2006]). The formation of a back in 
fast ICMEs is also supported by numerical simulations (see 
Figure 4 of Wu et al. [2005]). 

At the difference of previously studied MC-backs, the re- 
gion from '3' to '4' is compressed by the second MC. This 
region of interaction does not have the same characteristics 
as corotating interaction regions [CIRs, e.g. Pizzo, 1994] or 
typical merged interaction regions [MIRs, e.g. Burlaga and 
Ness, 1993, and references there in]. CIRs are due to a fast 
SW overtaking a slow SW, and they develop at larger dis- 
tance from the Sun (> 10 AU) than MIRs. The physics 
involved in the region between two interacting MCs is ex- 
pected to be different, since MCs have a moderate spatial 
extension and are structured by the magnetic field. It im- 
plies that the compressed plasma can be evacuated on the 
sides, and that the second MC is able to accelerate fully the 
first one. Without a significant reconnected flux between 
the two MCs (when the magnetic fields are nearly parallel) 
and after a transient period of time from the interaction, the 
two MCs are expected to travel together (a situation fully 
different than in the case of a fast SW overtaking a slow 
SW). 

The next substructure that can be identified is '4-7'. It 
is a very huge region (~ 0.3 AU) that presents a low vari- 
ance of B and clear signatures of an expanding MC, with 
a very strong magnetic field that rotates coherently (from 
north-east to south-west). This region shows an almost lin- 
ear decay of B (from B ~ 60 to B ~ 15 nT) consistent with 
the observed expansion, with a linear V profile from 900 to 
700 km s~^ during ~ 14 hours (equivalent to ~ 0.3 AU), 
a typical expansion rate observed in MCs [Demoulin et al., 
2008] . It also presents values of Tp significantly lower than 
Tex, low proton density (rip ~ 3 cm~^), and very low values 
of /3p (~ 10~^). All these signatures discard other well- 
known non-MCs structures (as e.g. corotating interaction 
regions) . 

Inside region '4-7', we also remark '5' and '6'. At '5' a 
strong decrease of Up starts and there are small variations 
in the large scale trend of 9b, 4>b, B, and V. At '6', /3p, Up, 
and Tp have peaks (here Tp reaches Tex). Thus, it is a priori 
unclear where to set the end of the second flux rope, we can 
distinguish three possible ends: '5', '6', or '7'. Finally, we 
set it at '7' because, at this position, a very strong MD is 
found and the coherence of B is lost. The presence of two 
flux ropes inside this ICME is also supported by applying 
the Grad-Shafranov technique [Hu and Sonnerup, 2002] to 
the plasma and magnetic field data at LI (Mostl, C, 2007, 
private communication) . 

An expanding structure with ICME signatures (e.g., low 
proton temperature, low pp, decreasing V and B profiles) 



still remains after '7'. At position '8', there is another MD 
(observed mainly in and (^b) together with a change in 
the decay rate of B. This MD is associated with a sud- 
den increase of Tp (and consequently of /3p, because rip re- 
mains approximately constant), while the velocity profile 
still presents roughly the same slope (i.e. the same expan- 
sion rate). Even when the expansion signatures (in both, 
V and B) and the coherence of B end at '9', low /?p and 
Tp < Tex remain along one hour, until '10'. Here B recovers 
values as low as those found in the typical solar wind, Tp re- 
covers the expected values (Tp ~ Tex), the coherence in the 
rotation of B is lost, and the level of fluctuations (e.g., in B, 
see panels 9b and 4>b) start to increase significantly. This 
large region, from '7' to '10', shows the characteristics of an 
extended back feature belonging to the second flux rope. 

The last panel of Figure 1 shows the geomagnetic effect 
of the ICME, as monitored by the Dst index (preliminary 
index downloaded from the OMNI data base). From '2' the 
intensity of the storm increases, because of the high intensity 
of the driven electric field {VB^ ~ 5 x 10* nT km s^^, with 
Bs being the southern component of interplanetary mag- 
netic field). 

The storm starts its recovery phase at ~ 08:00 UT, 
when the magnetospheric response is the strongest (reaching 
Dst = —263 nT); we notice a significant change in the decay 
rate after '3'. Later, Dst remains < —100 nT during part 

of the recovery phase (until '5'). Assuming a pure decay, 
decay times of r ~ 5 hours and r ~ 17 hours are found for 
the first (from the peak of the storm to '3') and second stage 
(beyond '3'), respectively. These two values are beyond the 
lower and close to the upper border of the typical range ob- 
tained for the ring current decay time, t = 14 ± 4 hours, 
by Dasso et al. [2002]. Larger decay times for the recovery 
phase have been associated with the presence of multiple IP 
structures near Earth [Xie et al., 2006]. 

2.2. First Magnetic Cloud 

In this section, we model the magnetic structure follow- 
ing the sheath (i.e., the data between '2' and '3') that we 
call the first cloud (MCI). 
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Figure 2. ACE observations (dots) of the magnetic 
field vector components in the cloud frame (see Sec- 
tion 2.2) for the first flux rope (observations between '2' 
and '3'). From upper to lower panels: -Ba,,cioud, By, cloud, 
and Bz, cloud. The solid lines correspond to the fitted 
curves obtained using the Lundquist's model. 
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To better understand the MC properties, we define a co- 
ordinate system linked to the cloud in which Scioud is ori- 
ented along the cloud axis (with -63, cloud > at the cloud 
axis). Since the speed direction of a cloud is mainly aligned 
with the Sun-Earth direction and is much larger than the 
spacecraft speed, we assume a rectilinear spacecraft trajec- 
tory in the cloud frame. This trajectory defines a direction 
d (pointing toward the Sun). Then, we define i/cioud in the 
direction Scioud x d and finally fcioud completes the right- 
handed orthonormal base (idoud, 2)cioud, Zcioud). 

We define the axis latitude angle (9) as the angle be- 
tween the cloud axis and the ecliptic plane, and the axis 
longitude angle [ip) as the one between the projection of the 
cloud axis on the ecliptic plane and the Earth-Sun direction 
(xose) measured counterclockwise [see Dasso et ai, 2006]. 
We also define the impact parameter, p, as the minimum 
distance from the spacecraft to the cloud axis. 

The local coordinate system is especially useful when p 
is small compared to the MC radius (R). In particular, for 
p = and a cloud described by a cylindrical magnetic con- 
figuration B{r) = Bz{r)z + B^{r)<j), we have ^doud ~ f and 
ydoud = 4> when the spacecraft leaves the cloud. In this par- 
ticular case, the magnetic field data will show: i3i:,cioud = 0, 
a large and coherent variation of By, cloud (with a change of 
sign), and an intermediate and coherent variation of Bz. doud, 
from low values at one cloud edge, achieving the maximum 
value at its axis, and returning to low values at the other 
edge. 

The minimum variance (MV) method [Sonnerup and 
Cahill, 1967] has been used to estimate the orientation 
of MCs [see e.g., Bothmer and Schwenn 1998, Lepping 
et al. 1990, Farrugia C. J. et al. 1999, Dasso et al. 2003, 
Guhsano et al. 2005]. In particular, using a synthetic set of 
ideal cylindrical clouds, Gulisano et al. [2007] have shown 
that the application of the MV technique to the normal- 
ized observed time series of the magnetic field {B{t) /\B{t)\) 
can provide very good estimations of the cloud axis (when 
< p/R < 0.7). This is feasible because when B{t) is nor- 
malized, the information of the rotation of the field vector is 
not mixed with possible changes in its absolute value during 
the observations. 

We apply the normalized MV method to the observations 
between '2' and '3'. We find a left handed flux rope, ori- 
ented such that 9mv~-15° and (y3MV = 125°, consistent with 
the SEN orientation determined in Section 2.1. From this 
orientation and the mean bulk speed (894 km s~^ for this 
range), we estimate the cloud size perpendicular to its axis 
as ~ 0.08 AU. 

In the MC frame (see Figure 2), the magnetic field has 
the typical shape observed in clouds. The almost con- 
stant -Bo;. cloud profile (with a mean value < -Ba:,cioud >~ 
— 11.5 nT), indicates that p is not zero. Moreover, a rough 
estimation of p can achieved using the expression found b y 
Guhsano et al. [2007] {p/R ~ y/{< Bo^.cioud > /Bo)/1.6), 
which gives p ~ (0.3 - 0.4)i?. 

The MC borders set as '2' and '3' are in agreement with 
the expected conservation of magnetic fiux across a plane 
perpendicular to Xdoud [Dasso et al, 2006], i.e., similar ar- 
eas below and above the curve Sy^cioud which equals before 
and after the MC center. These boundaries are also in agree- 
ment with the expected MD (current sheets) at the inter- 
faces between two structures with different connectivity (as 
the boundaries of a flux rope forming a MC). 

The MC axis, with lpmv = 125°, is dominantly pointing 
towards ycsE with a signiflcant contribution toward -xgse- 
This is consistent with the spacecraft passing through the 
right (west) leg of the flux rope and, thus, with the MC 
apex located toward the left of LI (observed from Earth 
to the Sun, with north upward). From the sign of the ob- 
served < -Ba:, cloud >, the cloud axis is toward the south of 
the ecliptic. 

Fixing the orientation provided by the MV method, we fit 
the free parameters of Lundquist's model (see Section 1 for 



the meaning of the free parameters) and obtain Bo = 57 nT 
and a = —40 AU~^. Solid lines in panels Bj,, cloud and 
Bz, cloud in Figure 2 show the curves obtained from the fit- 
ting, which are in a very good agreement with the observa- 
tions (dots). 

To validate our previous results, we perform a simulta- 
neous fitting (SF) of the geometrical and physical free pa- 
rameters using the same procedure and numerical code as 
Dasso et al. [2006]. We obtain 6Isf=-12°, yJsF=129°, Rsf = 
0.04 AU, p/RsF = 0.2, Bo.sF = 59 nT, asF = -42 AU'^ 
The agreement between the results of the SF and those from 
the MV method followed by the fitting of the physical cloud 
parameters is well within the precision of the methods (the 
differences are only ~ 4° in the orientation, 5% in the radius 
and 3% in Bo). Due to the radial propagation of the solar 
wind and due to a possible meridional stratification of the so- 
lar wind properties, distortions from a cylindrical cross sec- 
tion are expected [see, e.g., Liu et al, 2006b]. From a mul- 
tispacecraft analysis of the first MC observed by STEREO, 
an oblated transverse size with the major axis perpendicu- 
lar to the Sun-Earth direction was found [Liu et al, 2008b]. 
The good match between the observations and the cylindri- 
cal model (Figure 2) indicates that this flux rope could have 
an almost cylindrical configuration (a significant ellipticity 
would have consequences on the magnetic field profile, see 
Vandas and Romashets [2003]). Therefore, we conclude that 
the main characteristics of MCI are well determined. 

2.3. Second Magnetic Cloud 

In this section we model the second MC (MC2) observed 
from '4' to '7'. A decreasing velocity profile (see Section 2.1) 
indicates that the cloud is in expansion. The sudden change 
of Ob at '7' is the most significant of all changes observed 
in the field after '4'; furthermore, after this discontinuity, 
the expected coherence of B is lost. Thus, as previously dis- 
cussed in Section 2.1, we choose '7' as the rear boundary of 
MC2. 

The use of the MV method for this second MC does not 
provide meaningful results, since the impact parameter is 
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Figure 3. Magnetic field vector in GSE components 
from ACE data (dots) for the second MC (between '4' 
and '7') and fitted curve (solid line) using the model of 
Hidalgo et al. [2002b]. 
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very large as indicated by the low rotation of the magnetic 
field vector (Figure 1) and also because -B^:. cloud is the largest 
field component in the cloud frame [not shown]. Moreover, 
there is significant expansion, thus normalizing the field is 
not enough to fully remove this effect. The SF to Lundquist 
solution, even with a normalized field, cannot provide in this 
case a reliable result for the same reasons. When the impact 
parameter is low and the boundaries of the flux rope are well 
determined, flux rope modelling generally provides a good 
representation of the magnetic field configuration of a MC 
[Riley P. et al.. 2004]. liowever, models need to be tested us- 
ing simultaneous observations of different parts of the same 
fiux rope, e.g. as recently done using STEREO observations 
[Liu et al., 2008b] confirming the fiux rope geometry of the 
studied event. Since simultaneous observations from space- 
craft with a significant separation are not available for this 
event, we apply a different model and method for MC2, de- 
scribed below. This model has more freedom in comparison 
with Lundquist's model used before, such as the oblateness 
of the cloud cross section [Liu et al., 2006b]. 

Observations of some expanding MCs traveling in the 
solar wind are consistent with cylindrical expansion [e.g., 
Nakwacki et al, 2005; Dasso et al., 2007; Nakwacki et al., 
2008]. However, we can anticipate distortions from the cylin- 
drical shape for MC2 because of the presence of MCI in its 
front. Therefore, we compare observations with the model 
by Hidalgo et al. [2002b], which considers a magnetic field 
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Figure 4. GOES lightcurve from May 13 to 14. The 
fiares related to the eruptions we consider to be the 
sources of the two MCs (see Sections 2.2 and 2.3) are 
pointed by arrows at peak intensity (13:04 UT and 
16:57 UT for the first and second flare, respectively). 



with an elliptical cross section. This model has eight free 
parameters: the magnetic field strength at the cloud axis, 
the latitude (9) and longitude (tp) of the cloud axis, the 
impact parameter (p), the orientation of the elliptical cross 
section relative to the spacecraft path (i^), a parameter re- 
lated to the eccentricity of the cross section (77), and two 
other parameters related to the plasma current density. 

The fitted curves using the model of Hidalgo et al. [2002b] 
are in very good agreement with observations (Figure 3). 
The model corresponds to a left handed magnetic cloud with 
= (55 ± 10)° and (p = (120 ± 10)°. 

3. Solar clues for two source events 

3.1. The data 

In this section, we analyze solar data from the photo- 
sphere to the upper corona in search of two candidate solar 
events that could be the sources of the two MCs described 
in previous sections. 

The photospheric magnetic field evolution is analyzed us- 
ing observations from the Michelson Doppler Imager [MDI 
Scherrer and et al., 1995], on board the Solar and Helio- 
spheric Observatory (SOHO), which measures the line of 
sight magnetic field at the photosphere. These data are the 
average of 5 magnetograms with a cadence of 30 seconds. 
They are constructed once every 96 minutes. The error in 
the fiux densities per pixel in the averaged magnetograms is 
~ 9 G, and each pixel has a mean area of 1.96 Mm^. 

At the chromospheric level we have used full-disk Hq data 
from Big Bear Solar Observatory (BBSO) and the Optical 
Solar Patrol Network (OSPAN) at the National Solar Ob- 
servatory in Sacramento Peak. 

To identify changes in the low corona associated to 
the source regions of the identified MCs, the Extreme- 
Ultraviolet Imaging Telescope [FIT, Delaboudiniere J. -P. et 
al, 1995] on board SOHO is used. EIT images chromo- 
spheric and coronal material through four filters. In par- 
ticular, we have analyzed the 195 A band which images 
plasma at 1.5x10'' K. When available we have used observa- 
tions taken by the Transition Region and Coronal Explorer 
[TRACE, Handy, 1999] in the 171 A band. 

The identification of CMEs is done using LASCO on 
board SOHO, plus proxies for eruptions in the chromosphere 
and lower corona. For the analyzed time interval, LASCO 
imaged the solar corona from ~ 2 - 30 solar radii with two 
different coronagraphs: LASCO C2 (2 - 6 solar radii) and 
C3 (4 - 30 solar radii). 

3.2. Solar activity from May 11 to May 12, 2005 

Without any doubt, the main contribution to the larger 
cloud (MC2) and other IP signatures discussed in previ- 
ous sections is the halo CME appearing in LASCO C2 at 
17:22 UT on May 13 (see Sections 3.4 and 3.5). However, 
since we have found that the ACE magnetic field and plasma 
observations from 05:42 UT on May 15 ('3') to 10:30 UT on 
May 17 ('10') can be interpreted as being comprised by two 
different structures (see Sections 2.2- 2.3), we describe the 
solar activity observed by EIT and LASCO preceding the 
major CME on May 13. 

From May 11 until May 13, 2005, solar activity is mainly 
concentrated in two ARs, AR 10758 located in the south- 
ern hemisphere and AR 10759 in the northern hemisphere, 
which produces the most intense events. Activity progres- 
sively increases in the later region from May 11 to 12; then, 
fiares reach level 2B in Ha and Ml. 6 in soft X-rays. 

On May 11 a CME is first seen in LASCO C2 at 20:13 UT 
above the SW limb. This event is classifled as a full halo. 



DASSO ET AL.: TRACKING TWO CONSECUTIVE MAGNETIC CLOUDS 



X - 7 



2005/05/13 12:48:03 UT 



i 



2005/05/13 12:54:02 UT 



Figure 5. OSPAN Ha images of the flare at 12:49 UT. The left panel shows the portion of the filament 
lying along the E-W northern magnetic inversion line expanding northward. Flare ribbons and flare loops 
linking the ribbons are visible in the right panel. The E-W middle portion of the filament is no longer 
visible. 
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Figure 6. Evidence of flux cancellation at the northern E-W portion of the filament channel in AR 10759. 
Left: BBSO Ha image early on May 13 showing the non-uniform shape of the E-W filament portion. The 
two thick white arrows indicate the axis orientations of the two identified MCs projected on the solar 
disk. Center: MDI magnetogram closest in time to the flare at 12:49 UT (the field has been saturated 
above/below ± 100 G). The arrow points to a small positive flux concentration that moves towards the 
negative held across the filament channel. Several of such fiux cancellation episodes occur during May 13. 
Right: An overlay of a flare image and the ± 50 G MDI isocontour (white positive field, black negative 
field). In all images the axes are in arc sec. 



The linear and second order fittings to LASCO C2 ob- 
servations give plane-of-sky (POS) speeds of 550 km s~^ 
and 495 km s~\ respectively (from LASCO CME Cata- 
log, http://cdaw.gsfc. nasa.gov. /CMEJist/index.html). An 
Ml.l flare in NOAA AR 10758, located at S11W51, start- 
ing at 19:22 UT can be associated with this event. Sub- 
frames from the EIT shutterless campaign show a raising 
feature in absorption by 19:22 UT. However, the closeness 
in time/space between the two MCs discussed in Section 2.1 
makes it highly improbable for this event to be the source of 
MCI, since by simple assumption of constant speed it would 



take ~ 3.1 days to reach Earth and then would be expected 
to arrive before May 15, i.e. earlier than MCI. 

On May 12, an event (flare, dimming, and loops in expan- 
sion) is seen in EIT at about 1:57 UT. This event originates 
in AR 10759 and is accompanied with a C9.4 flare located 
at N11E31 with peak at 1:10 UT. No CME is observed in 
LASCO. However, a diffuse and gusty outward flow is ev- 
ident in LASCO C2. There is another apparent event in 
EIT at about 13:56 UT, to the south of AR 10760. The re- 
lated C3.0 flare occurs in S11W62 at 13:40 UT. This event 
shows an association with a CME, weak and slowly travel- 
ing in the SW direction. This CME has a strong component 
in the POS and, therefore, we do not expect it to be sig- 
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nificantly directed towards Earth. Furthermore, from the 
orientation of AR 10760 polarities during their disk transit, 
we can infer from MDI data that its magnetic helicity sign 
is positive, contrary to that of both studied MCs [see the 
description and interpretation of magnetic tongues in Lopez 
Fuentes et al., 2000, 2003]. 

After discarding these previous events, and considering 
the closeness in time/space between the two MCs analyzed 
in Section 2, we concentrate in the activity observed during 
May 13. In the next sections we will present observations 
and discuss, in particular, two flares and associated filament 
eruptions that we consider to be the solar source events of 
the clouds. The flare timings are given by GOES soft X-ray 
data (Figure 4). 

3.3. May 13 events: Photospheric and chromospheric 
signatures 

Along May 13 the filament in AR 10759 is seen to activate 
several times, mainly in its northern fraction. At 12:49 UT 
a CI. 5 flare starts in N17E15, with peak at 13:04 UT. The 
rising filament is seen in Ha and also in TRACE (see the 
left Figure 5 and Section 3.4). At the chromospheric level 
a two-ribbon flare, classified as a subflare faint, develops 
along the magnetic inversion line at the north of AR 10759 
(Figure 5). This is the only event showing eruptive char- 
acteristics at that time on the solar disk. This portion of 
the inversion line is oriented in the E-W direction. Part of 
the filament extending along this E-W inversion line is no 
longer visible by 12:54:02 UT. Flare loops linking the two 
ribbons are also visible in Hq at this time. What may be 
the cause of the dcstabilization of the E-W portion of the 
filament? We have analyzed the magnetic changes observed 
in the filament channel previous to the flare at 12:49 UT, for 
this we use MDI line of sight magnetic maps. Continuous 
flux cancellation is observed at the magnetic inversion line 
along which the erupting part of the filament lies. The E-W 
northern fraction of the filament shows a non-uniform shape 
(as being formed by several sections) and some barbs or feet 
are present (Figure 6 left). A small positive polarity is seen 
to intrude in the negative field in between two sections of 
the filament (sec the arrow in Figure 6 center); this single 
intrusion implies a cancellation of 5 x 10^'^ Mx. The close 
relationship between the fiarc and the magnetic cancellation 
site is shown in the right panel of Figure 6. MDI low tempo- 
ral and spatial resolution does not allow us to follow in detail 
the evolution of the small flux concentrations and ascertain 
the clear association between their cancellation and filament 
eruption; however, we believe this to be the most plausible 
cause. Small-scale magnetic changes in flux concentrations 
along filament channels have been frequently reported a few 
hours prior to local filament restructuring [sec e.g., Chae 
et al, 2000; Deng et al, 2002; Wood and Martens, 2003; 
Schmieder et al, 2006]. Furthermore, Zhang et al [2001] 
attributed the origin of the major solar flare and fllament 
eruption on July 14, 2000, to flux cancellation at many sites 
in the vicinity of the active region fllament. 

Taking into account that the CI. 5 flare occurs by only 
about 4 hours earlier than the main flare and eruption in 
AR 10759 (see below) and the similar orientation of the 
MCI axis and of the magnetic inversion line along which 
the erupting filament lies (Figure 6 left), we consider that 
the partial eruption of the AR filament is the source of this 
first cloud. There is evidence for the breakup of filaments 
into more than one segment, [see e.g. Martin and Ramsey, 
1972] and a recent example in Maltagliati et al [2006]. Even 
some filaments erupt only partially, or saying it in a differ- 
ent way different portions of the same filament may erupt at 
different times and trigger different fiares [Martin and Ram- 
sey, 1972; Tang, 1986; Maltagliati et al, 2006; Gibson et al. 



2006a; Gihson and Fan, 2006; Liu et al, 2008a]. This kind 
of eruptions have been discussed in the frame of the eruptive 
fiux rope model for CMEs [see the review by Gibson et al, 
2006b]. Moreover, as discussed below in Section 3.5, on May 
13 there are no other evidences of CMEs in LASCO until 
17:22 UT. 

The solar event following the one just described is the 
major M8.0/2B flare that starts at 16:13 UT and peaks at 
16:57 UT. This flare ended on May 14 at ~ 17:00 UT, be- 
ing a long duration event (see Figure 4). The photospheric 
and chromospheric observations corresponding to this event 
have been shown and discussed by Yurchyshyn et al [2006] 
and Liu et al [2007]. The former authors have associated 
this flare and eruption to a single MC comprising the two 
we have identified in this paper. We want to stress here that 
the direction of MC2 axis is in agreement with the orienta- 
tion of the N-S fraction of the AR filament that lies along 
the main magnetic inversion line (see Figure 6, left panel). 

3.4. May 13 events: The lower corona signatures 

In the lower corona, TRACE observes the event at 
12:49 UT in 171 A (Figure 7). Images before 12:39:49 UT 
and after 13:04:27 UT are not useful because of a strong 
particle shower. TRACE data are dominated by a large- 
scale sigmoidal structure (with an inverse S shape) previous 
to the flare weU visible by ~ 12:40 UT. The small scale 
loops linking the two Ha ribbons are not distinguishable in 
TRACE images. By 12:49:43 UT the filament can be ob- 
served as a dark curved structure surrounded by a brighter 
one, probably formed by heated filament material (see Fig- 
ure 7). This bright structure is seen expanding upward until 
13:04:27 UT, though it is hard to assure that it erupts be- 
cause of lack of visibility in the following images (Figure 7, 
central panels). 

Later on, by ~ 16:00 UT (last panel in Figure 7) and be- 
fore the M8.0 fiare, TRACE images appear again dominated 
by the large-scale sigmoid. This structure is similar to the 
one described by Liu et al [2007] seen before the M8.0 flare; 
in their Figure 1 they show the diflferent parts that consti- 
tute the sigmoid: the magnetic elbows, the envelope loops 
[according to the nomenclature of Moore et al, 2001] and 
the footpoints of the sigmoidal fields. Thus, the rising of 
the bright curved structure, shown in the central panels of 
Figure 7, is most likely linked to the filament eruption during 
the first and less intense two-ribbon flare that we associate 
with MCI. Yurchyshyn et al [2006] have suggested that the 
bright structure seen in the central panels is the core of the 
sigmoid that later starts a slow ascent (by ~12:55 UT) and 
become active during the M8.0 flare, but in fact it erupted 
earlier and reactivated later. 

The major eruptive event, linked to an M8.0 flare, is visi- 
ble in EIT images starting from 16:37 UT (see Figure 8). It 
exhibits several CME-associated phenomena [Hudson and 
Oliver, 2001]: coronal dimmings, a large-scale EIT wave 
and a post-eruption arcade. As it was pointed out by Liu 
et al. [2007] , the darkest part of the dimming comprises two 
patches to the NE and SW of the erupting active region. 
The SW dimming is associated with a transient coronal hole 
(shown by an arrow in Figure 8b). It was suggested that 
transient coronal holes can be interpreted as footpoints of 
the ejected flux rope [sec, e.g., Webb et al, 2000; Mandrini 
et al, 2005]. The absence of the clearly visible second tran- 
sient coronal hole in this event may indicate that the second 
leg of the erupted flux rope may have already reconnected 
with the ambient coronal magnetic field [this reconnection 
changes the very extended field lines rooted in the dimming 
by closed magnetic connections Attrill et al, 2006; Kahler 
and Hudson, 2001]. The overall morphology of twin dim- 
mings and the presence of the post-eruption arcade (Fig- 
ure 8b) suggest that an eruption of a coronal magnetic flux 
rope had taken place [e.g., Hudson and Oliver, 2001]. 
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Figure 7. TRACE images in 171 A. The first tliree panels sliow different stages of the expansion of 
tlie filament embedded in a brighter curved structure during the flare at 12:49 UT. The last panel shows 
that the coronal loop structure recovers its global shape after the CI. 5 flare and before the M8.0 flare at 
16:32 UT. The same region is showed in the four panels, though the TRACE field of view is shifted for 
the image at 16:08:43 UT. 




Figure 8. (a) SOHO/EIT difference image in the Fe XII bandpass (195 A) showing a large coronal dim- 
ming (dark areas) corresponding to the CME associated with the M8.0 flare on May 13, 2005. The last 
pre-emption image taken at 16:26 UT was subtracted from the image taken at 17:07 UT. (b) SOHO/EIT 
image in the Fe XII bandpass (195 A) taken at 22:57 UT on May 13, 2005 (after the start of the CME 
associated with the M8.0 flare). Post-eruption arcade and a transient coronal hole are marked with top 
and bottom arrows respectively. 





Figure 9. LASCO images of the halo CME on May 13, 2005, as recorded by a) LASCO C2 and b) 
LASCO C3. A prc-cvcnt image has been subtracted in both cases. These are the two only full-field of 
view images taken of the event. The arrows indicate the edge of the faintest front. 



3.5. May 13 events: The upper corona observations 

On May 13, 2005, a full halo CME with brightness asym- 
metry towards the NE first appears in the field of view of 
LASCO C2 at 17:22 UT (see Figure 9a). Five more partial 
snapshots of the halo CME were captured by LASCO C3 
from 17:20 to 17:50 UT. The closeness in time of the lat- 



ter images does not allow to see major changes from one to 

another, so we present here only one of them (Figure 9b). 
The estimated POS speed for the fastest front, as calculated 
by the LASCO CME catalog, is ~1690 km s"^ Further in- 
formation on the kinematics of this CME and its associated 
shock can be found in Section 4 and in Reiner et al. [2007] . 
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Figure 10. Dynamic spectral plot showing the radio emission recorded by the WAVES/TNR receiver 
on board Wind between May 13 at 12:00 UT and May 15 at 12:00 UT. The vertical axis represents 1// 
[kHz~^] while the horizontal axis corresponds to time. Three identified type II emissions are labeled (the 
entity (a) is associated to (1) and (2) which arc the fundamental and harmonic emission at the plasma 
frequency, and the entity (b) is associated to (3)). The black arrow in the horizontal axis denotes the 
shock arrival time. Tick marks on the right vertical axis indicate the location of TNR frequency channels 
for the plotted frequency range. 



Enhancement of the LASCO images presented in Figure 9 
allows to distinguish a very faint edge toward the south and 
north, about a couple of solar radii ahead of the bright- 
est one (indicated by arrows). This feature suggests the 
presence of a shock, as described by Vourlidas et al. [2003] 
and Vourlidas [2006] . Finally, the CME related to the solar 
eruption at 12:49 UT is probably so faint, that its detection 
was not likely due to sensitivity limitations of LASCO [e.g., 
Tripathi et al,, 2004; Yashiro et al., 2005]. Furthermore, 
Yashiro et al. [2005] conclude that all fast and wide CMEs 
are detectable by LASCO, but slow and narrow CMEs may 
not be visible when they originate from disk center. In ad- 
dition to this restriction, which impedes the visualization of 
a CME if it is feeble, the LASCO cadence during and after 
the CI. 5 flare at 12:49 UT was of 30 minutes, as opposite to 
the usual 12 minutes of LASCO C2. 

4. Two structures in the inner heliosphere: 
Remote radio type II observations 

When a halo CME leaves the field of view of a corona- 
graph, it is occasionally possible to further track its evo- 
lution by means of measurements in the radio regime of 
the electromagnetic spectrum. The shock commonly driven 
by a CME excites electromagnetic waves emitted at the lo- 
cal plasma frequency (/ ~ and/or its harmonic. As 
the shock travels outward from the solar corona, the local 
plasma frcqucncj' decreases, and a drifting signal, a type II 
radio burst, is generated. This brand of emissions can be 
followed drifting in frequency from ~ 150 MHz down to 
^ 25 kHz, which is the approximate value of the local plasma 
frequency in the vicinity of Earth [Leblanc et al., 1998]. 
Then, the shock traveling toward Earth is detected by in- 
struments in the IP space, such as WAVES on board Wind 
[Bougeret J. L. et al, 1995]. The ICME can also be de- 
tected by ground-based antennas measuring interplanetary 
scintillation of radio sources [e.g., Manoharan, 2006]. 

The metric components of type II bursts have long been 
used by space weather forecasters as input to models of 



shock arrival time prediction [e.g.. Dryer and Smart, 1984; 
Smith and Dryer, 1995; Fry et al, 2001]. Nevertheless, there 
is ongoing debate on the mechanism that generates them, 
questioning their direct relationship with shocks ahead of 
CMEs [e.g., Pick, 1999; Gopalswamy et al, 2001a; Clafien 
and Aurass, 2002; Cane and Erickson, 2005]. Therefore, 
Cremades et al. [2007] attempted to improve space weather 
forecasts by employing information of type II radio bursts 
in the kilometric domain (30 — 300 kHz), which are indeed 
believed to be generated by IP shocks [Cane et al, 1987]. 
Accordingly, the slope of the drifting radio emission can be 
taken as a proxy for the associated shock speed [Reiner 
et al, 1998]. Furthermore, for each point of emission de- 
tected in the space-borne radio receivers, it is possible to 
obtain a proxy of the radial distance from the Sun through 
density models [e.g., Saito et al., 1977; Leblanc et al., 1998; 
Hoang et al., 2007]. The technique used to characterize the 
temporally-associated radio signals makes use of the density 
model developed by Leblanc et al. [1998]. To obtain the ap- 
proximate sites of the emissions, and consequently a proxy 
of the shock speed, it is only needed to enter into the density 
model an estimate of the plasma density near Earth, usu- 
ally between 5 and 9 cm^"^. The top frequency taken into 
account as input for the technique (300 kHz for the funda- 
mental emission) corresponds approximately to a distance of 
20 solar radii according to that heliospheric density model. 

We assume that at those heights most of the decelera- 
tion has already taken place, and thus that the slope of the 
linear profile in Figure 10 is enough to characterize the kine- 
matics of a shock that travels at a constant speed. Beyond 
20 solar radii, the corresponding frequency channels of the 
WAVES/RADl detector arc scarce in comparison with those 
of WAVES/TNR (see rationale at the end of this Section). 
That is why the main goal of the latter detector is the track- 
ing of radio phenomena up to the orbit of Earth. Further 
details on the technique, its application and success relative 
to other methods can be found in Cremades et al. [2007]. 

After the solar events on May 13, a number of type II 
radio features in the kilometric domain were observed by 
the Thermal Noise Receiver (TNR) experiment on WAVES. 
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Figure 10 displays the radio intensity between May 13 at 
12:00 UT and May 15 at 12:00 UT, and from / -250 kHz to 
/=20 kHz. In addition, each of the identified low-frequency 
type II emissions has been labeled as follows: (1) a radio 
emission drifting at a speed of 1400 km s~^, (2) an emission 
drifting at a speed of ~1450 km s~^ and certainly being the 
harmonic of the first feature, and (3) a slower radio emis- 
sion drifting at a speed of ~850 km s~^. The application of 
the technique described in Cremades et al. [2007] to emis- 
sion (1) yielded for its associated shock a predicted arrival 
time at ~02:40 UT on May 15, which is in good agreement 
with the real shock arrival at Earth (02:11 UT). The arrival 
time derived from radio signal (2) yielded ~23:20 UT on 
May 14. The error of ~3 hours might be attributed to the 
'patchiness' of the emission. This effect arises because a sin- 
gle shock may drive type II bursts at multiple places along 
the shock surface, with the possibility of a large non-radial 
component of the shock velocity. In addition, there could be 
significant differences in the drift rates if the density gradi- 
ent or the local shock speed varies along the shock surface. 
Under these circumstances entities (1, 2) and (3) could be 
the manifestation of the shock driven by only one ICME. 
However, we have no data to justify any of these hypoth- 
esis. Moreover, the large difference in propagation speed 
between entities (1, 2) and (3) suggests that they must be 
indeed produced by two distinct structures traveling in the 
solar wind or by the same structure traveling with a sig- 
nificantly slower speed at a first stage of its journey and 
with a faster speed at its second stage. We will interpret 
this second possibility in the frame of the global observa- 
tions in next section (Section 5). Finally, the prediction of 
the arrival time of entity (3) yielded ~10:20 UT on May 15. 
The linear backtracing of this last emission in time matches 
well with the earlier solar event on May 13. On the other 
hand, entities (1) and (2) show good temporal association 
with the occurrence of the second solar event on May 13 (at 
16:57 UT). 

Based on the observations of WAVES/TNR we present 
here a different interpretation of the sources of the radio 
signatures than Reiner et al. [2007] , who used radio data ac- 
quired by the WAVES/RADl detector (see next paragraph). 
The TNR observations indicate the existence of 2 interplan- 
etary shock entities, which from now on will be referred to 
as entity (a) -the shock traveling at ~1400 km , corre- 
sponding to (1) and (2)- and entity (b) -the shock traveling 
at 850 km s~^, related to (3). Entity (b) presumably draws 
near entity (a) around 07:00 UT on May 14, at an approxi- 
mate distance of 98 solar radii from the Sun. There are no 
visible signs of 'cannibalism' [Gopalswarny et al., 2001b] in 
radio data, likely due to the weakness of entity (b) and/or 
because there is no significant reconnection between the two 
entities. Moreover, it is worth to note that entity (a) was 
not visible before 06:00 UT on May 14, as opposite to the 
slower entity (b). A plausible reason for this effect is that 
entity (a) may have been moving in an environment plasma 
faster than the typical solar wind, therefore without leading 
shock formation. 

The kinematics analysis performed for the radio event on 
May 13-15, 2005 by Reiner et al. [2007] takes into account 
entity (a) only. Those authors present spectral plots of the 
event that reach down in frequency up to the coverage of the 
WAVES/RADl detector. As seen in Figure 10 of that study, 
entity (b) is not obvious from RADl measurements. How- 
ever, it does show up in TNR data, suggesting that Reiner 
et al. [2007] did not consider those data in their analysis. It 
is most likely that the poor frequency resolution of RADl 
below 256 kHz is guilty of overseeing entity (b). The RADl 
detector employs at a time only 32 frequency channels to 
monitor the range 20-1040 kHz, using only 21 channels be- 
low 256 kHz and interpolating to obtain the rest of the data. 
This may lead to false apparent broadband emissions and to 



the oversight of some narrow-band features in the RADl do- 
main. Conversely, the TNR receiver employs 96 channels to 
cover the frequency range of 4-256 kHz, over five logarith- 
mically spaced frequency bands [Bougeret J. L. et al., 1995, 
also see location of frequency channels as tick marks on the 
right vertical axis of Figure 10]. This configuration achieves 
better frequency resolution at distances greater than ~20 
solar radii, and hence makes it best suitable to study emis- 
sions beyond those distances. 

5. Conclusion 

We propose here a fully consistent physical scenario for 
the chain of events from May 13 to 17, 2005, which is based 
on a detailed analysis of the observations presented in pre- 
vious sections. 

Two solar events occurred on May 13 with a difference 
of ~ 4 hours, both from AR 10759. The first one, a CI. 5 
two-ribbon flare that was associated with the ejection of a 
part of the AR filament extending along the E-W portion of 
the inversion line (north of the AR) . The second and largest 
one, an M8.0 two-ribbon associated with the ejection of the 
filament lying along the N-S portion of the inversion line 
and a CME observed in LASCO C2. 

In the IP medium we identify an anomalously long ICME 
(more than 2 days). Following the results of Burlaga et al. 
[2001, 2003], this suggest the presence of two flux ropes (see 
Sect. 2.1). Indeed, we have found two clouds (MCI and 
MC2). They are stacked together and only the latitude, 6b, 
of the magnetic field permits to identify in the data two 
coherent structures separated by a region (called '3-4' in 
Fig. 1) that have a different temporal evolution from neigh- 
boring regions. Such conclusion is supported by the various 
attempts done to fit the data by one or two flux ropes. In- 
deed, only the two flux rope models presented in Section 2 
give a reasonably good fit to the data. Figure 11 shows a 
sketch of the two clouds near Earth, in which we have taken 
into account their orientations, relative sizes, and the space- 
craft trajectory during the observations. The two MCs are 
oriented such that: 9mci = (^12±5)°, ipMci = (125±5)°, 
Omc2 = (55±10)°, and <fiMC2 = (120±10)°. MCI presents 
a low impact parameter {p/R ~ 0.2 — 0.3) and a radius of 
~ 0.04 AU, while MC2 is rrmch larger and presents a larger 
impact parameter; both MCs have negative magnetic helic- 
ity. MCI is small, has a very intense magnetic field, the 
signature of an overtaking flow penetrating deep in the out 
bound, and almost no portion with low proton temperature 
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Figure 11. Schematic global view of the two MCs and 
the trajectory of the spacecraft that observed them. The 
first MC is drawn with a continuous line (MCI), while 
the second one with a dashed line (MC2). 
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(sec Figure 1); all these observations can be interpreted as 
the consequence of a compression by MC2. 

The left panel in Figure 6 shows the projection on the 
solar surface of the axes of MCI and MC2 and the sense 
of their azimuthal field component, for comparison with the 
orientation of the erupting parts of the filament associated 
with the CI. 5 and M8.0 flares, respectively. Only the sec- 
ond flare has an observed associated CME, and it is apriori 
surprising to associate a CI. 5 flare to a fast ejection (with 
a velocity ~ 850 km s^^ from radio data); however, it has 
been shown by Gopalswamy et al. [2003, 2005] that the cor- 
relation between CME speed and flare intensity is weak. 
Apart these two issues, the other characteristics of the solar 
and interplanetary events are in agreement as follows. The 
erupting fllament orientations and directions of the field in 
their associated arcades (see right panel of Figure 5 and Fig- 
ure 8b) are in good agreement with the axial and azimuthal 
field components in the clouds. The AR field, as well as the 
two MCs have negative magnetic helicity. Altogether, with 
the transit time in the expected time interval, this demon- 
strates the associations between the two solar eruptions to 
the two interacting MCs detected near Earth. 

Moreover, an evidence of the interaction is directly found 
in the interplanetary radio data, as follow. Radio type II 
observations (Section 4, entity (b)) are consistent with a 
speed for the shock driven by the first ejecta (MCI) of 
~ 850 km s~^, when it travels from the Sun up to ~ 0.5 AU 
(from ~13:00 UT on May 13 to ~07:00 UT on May 14). The 
shock driven by the second ejecta (MC2) seems to be weaker 
and is not detected using radio data in this range of time, 
probably because MC2 travels in a non-typical solar wind 
medium, i.e., a medium perturbed by the previous passage 
of MCI, so having a faster speed. It is also possible that this 
second shock starts to interact with the trailing part of MCI 
at earlier times just after its eruption and, thus, it turns to 
be weaker, as shown from numerical simulations of interac- 
tion of fiux ropes [Xiong et al., 2007]. Then, at times when 
MC2 reaches MCI (at ~ 0.5 AU, on May 14 at ~07:00 UT), 
when the exchange of momentum between the two clouds is 
most efficient, MCI is accelerated and MC2 is decelerated. 
Then, both start to travel at similar speeds ~ 1400 km , 
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Figure 12. Cartoon representing the position (as dis- 
tance from the Sun along the Sun-Earth line) of the 
events ejected from the Sun in function of time. Verti- 
cal dotted lines (ticks at SI and S2) represent the time 
of beginning of solar fiares associated with both ejections 
(May 13, ~ 13:00 UT and ~ 16:00 UT, respectivelly. 
Vertical dashed fines (ticks at Shock, MCI, and MC2) 
represent the arrival of the shock (May 15, 02:11 UT) 
and of the center of magnetic clouds 1 and 2 (May 15, ~ 
09:00 UT and - 21:00 UT, respectively). 



but keeping their individuality. The shock wave driven by 
this new combined structure is the cause of the radio emis- 
sion observed as entity (a), identified in Section 4. During 
the second half of the transit to Earth, MCI is compressed 
from behind by MC2. 

Summarizing, we propose the following interpretation of 
the sequence of events that occurred on May 13, 2005: two 
ejective solar events occurred, which could be tracked to the 
Earth environment and were observed as two attached, but 
non-merged, magnetic clouds. Figure 12 shows the position 
versus time (along the Sun-Earth line) of both interplane- 
tary events, as inferred from the sequence of the observations 
used for our analyses from the Sun to 1 AU. The presence 
of structures preserving multiple fiux ropes, as the multi- 
ple structure studied here, has been observed in some cases 
[Wang et al., 2002]. Signatures of the interaction between 
ejecta, similar to those found in this study, have also been 
identified from observations with the following properties: 
(i) acceleration of the leading ejecta and deceleration of the 
trailing one, (ii) compressed field and plasma in the lead- 
ing ejecta, (iii) weak shock or disappearance of the shock 
driven by the second ejecta, and (iv) strengthening of the 
shock driven by the leading one [Farrugia and Berdichevsky, 
2004]. 

This study especially illustrates the need of combining so- 
lar, in situ and remote sensing of interplanetary propagation 
of the ejecta, in order to reveal the physical processes. Pre- 
vious studies of the same event, but mostly limited to one of 

these three domains have concluded on the presence of only 
one ejecta (with varying characteristics depending on the 
study). In the present case, the peculiar characteristics of 
the in situ observed magnetic field first alert us on the pos- 
sible presence of two MCs. We confirm this by performing a 
detailed modeling of the magnetic field exploring various flux 
rope boundaries. This stimulates a deeper study of data in 
the two other domains, with the new perspective of testing 
two ejecta. The solar data were crucial to first confirm the 
presence of two eruptions, and to define the timing and the 
spatial organization of the erupting magnetic configurations. 
This brought us to carry a deeper analysis of interplanetary 
radio observations. Indeed, two type II bursts, traveling 
at different velocities, were found. These observations also 
permit to localize where the interaction occurred (about half 
way between the Sun and the Earth). Then, all the pieces 
of the puzzle finally fit together, re-enforcing our interpreta- 
tion separately in each domain. However, the possibility of 
a single complex event (as proposed by [Yurchyshyn et al., 
2006; Reiner et al., 2007]) cannot be fully discarded and our 
proposed scenario, based on compelling evidence, must be 
further investigated. 

Such studies of interacting ejecta will greatly benefit from 
observations taken with spacecraft in quadrature, as already 
done for some events [e.g., Rust et al, 2005], and as available 
with STEREO spacecrafts. This has already be done with a 
spacecraft configuration far from optimum conditions [Har- 
rison R. A. et al., 2008]. On top of providing in situ mea- 
surements of physical parameters at three locations (com- 
bining ACE with STEREO) and monitoring type II bursts, 
there are imminent new possibilities of imaging ICMEs in 
their journey through interplanetary space. Then, we will 
have the opportunity to derive precise constraints on the 
physical mechanisms of interaction of an ejecta with a faster 
one, or simply with a fast overtaking solar wind. 
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